RU9552 



Exact methods for Campi plots 

K. C. Chase and A. Z. Mekjian 

Department of Physics, Rutgers University 
Piscataway, New Jersey 08855 
(February 9, 2008) 

Abstract 

We introduce for canonical fragmentation models an exact method for com- 
puting expectation values which exclude the largest cluster. This method 
allows for the computation of the reduced multiplicity and other quantities of 
interest introduced by Campi, and a comparison shows that the percolation 
model and a recent canonical model differ mostly only in small respects in 

these ensemble averages. 
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Campi and Krivine [ljjj] introduced a method for distinguishing fragmentation models 
from one another. By comparing various expectation values in which the largest cluster 
is excluded but the particle number and fragment multiplicity are held fixed, they showed 
that the percolation model has a distinctly different behavior than many competing nuclear 
fragmentation models. 

In this paper, we analyze another statistical weight we have been using recently and show 
that it shares many of the same properties as percolation theory, a point already apparent 
from a consideration of critical exponents f3|. The method used in this paper is unusual in 
that it is an exact computational method. Monte Carlo sampling is avoided by exploiting 
some properties of the partition functions, and as such there are no statistical errors, an 
immense improvement over earlier methods. 

We begin by assuming that each fragmentation outcome happens with a probability 
proportional to the Gibbs weight 



W(n) = n ^ , (1) 
k>l nk - 

where n k is the number of fragments of size (or charge) k, and x k is a parameter associated 
with k sized fragments. We then define the microcanonical partition function as 

Zf\x)= £ Win), (2) 

where 7r m (A) is the set of partitions of A nucleons into m fragments, i.e. J2k^n k — A, 
J2k n k = m - These partition functions satisfy the identity 

which allows for the computation of the partition functions recursively from J2k{ n k) — m , 
since 

(n k )= Xk £-lnZ^= Xk ^ r . (4) 

Campi defined reduced moments as moments in which the largest cluster is excluded 
from the measure, i.e. 



M s (n) = 5> s n fc - ^ (5) 
k 

where fc max is the size of the largest cluster. This suggested the definition of the reduced 
variance 72 |],|2] for a fragmentation event should be 

M 2 M gfc k 2 n k - fc max 

Its expectation value can be computed by breaking events into classes specified by fc max , and 
summing the expectation values over those classes with the appropriate weight, i.e. 

( 72 > = (m - 1) £ Pr(fc m ax) Efcfc2 | n ; ) _! / T aX \2~ A:Lx > ( 7 ) 

"■max N ' 

where Pr(fc max ) is the probability of fc max being the largest cluster size and (nfc)(fc max ) is the 
expectation value of n k when fc max is fixed. 

To compute expectation values in which the largest cluster size is fixed we need to 
compute the partition function for such ensembles. Clearly this partition function is given 
by all the terms in the microcanonical partition function which have the highest x k 

in the term. Consider 

-Z < f\x 1 ,...,x km ^ 1 ,0,...,0) . (8) 

We see that AZ^™' (^ max ) is the partition function for ensembles with fixed maximum cluster 
size /c max , since the first term collects all terms with £fc max or lower, and the second term 
eliminates those terms which don't have an £/c max - From this result we can determine Pr(fc max ) 
and (nfc)(/c max ), which are given by 

>M*-) = iff"' ^ , 0>) 



(w fc )(fc n 



k > k ma ,x 

gi"!fc 1) (xi,...,Xfc m „,0,-,0) , _ , nn x 

A7(M/L \ 



This method is quite general and can be applied to other models. For example, equiparti- 
tioning models, which have weights given by 

W(n) = J] *T (11) 
fc>i 

can also be analyzed by this method with some minor modifications. For example, x& = 1 
is the model used by Sobotka and Moretto |§. 

With these identities there is sufficient information to compute (72) and other reduced 
moments for any Xk- We use Xk = x/k T for a variety of reasons discussed elsewhere |§. 
Campi and Krivine ||] following Mekjian Q considered this model with r = 1.0 and showed 
that its reduced variance and other related expectation values had a distinctly different 
behavior than the percolation model. Plotting the expected reduced variance vs. (m — l)/A, 
(72) (m) has a single peak. The location, height and width of this peak for the two models 
(and other models they considered) are completely different, suggesting the usefulness of this 
plot in distinguishing fragmentation models. The choice r = was considered by Gross, 
et. al. |T0|-|r2|], and a different model was analyzed by Pan and Das Gupta |T3]. 



Since that time our interest has turned to the choice r = 2.5 because of similarities with 
percolation theory and Bose condensation. Namely, the sudden appearance of an infinite 
cluster in the infinite A limit and the presence of condensation phenomena. As such, we 
have recomputed the Campi plots for this model and have discovered that they duplicate 
the percolation model results in many respects. Figure |I|(a) shows the results. The height 
and location of the peak of the reduced variance (72) are the same in both models. The only 
significant difference is the width of the peak which is larger in the Gibbs model than in the 
percolation model. Plots of (fc max )(m) vs. (m — l)/A given in Fig. |](b) are also very similar 
for both models, and the scaling behavior of the position, width and height with changing 
A also agree. 

Another plot suggested by Campi [13] is to divide the event space by the maximum 



cluster size of the event A; max an d the reduced second moment M 2 and plot the probability 
of the canonical model being at any particular point on the graph. This can be also be done 



exactly for Gibbs models in a way completely analogous to the way given above. Define 
the partition function Z A (rri2',x) as the sum of the Gibbs weight Eq. ([I]) over all partition 
vectors n which satisfy J2k ^ n k — A Y^k k 2 n k = mi. This can be computed by the following 
recursion, 

Z A (m 2 ;x) = —^kx k Z A - k (m 2 - k 2 ;x) (12) 

k 

with Zo(rri2; x) = <5 m2j o- If we define AZ as in Eq. we find again the partition function 
conditioned on /c max being fixed, which is proportional to the probability of having an event 
with both m2 and k max . Figure [| plots this probability profile as a contour plot, which 
reveals the events are centered on a particular region in this phase space. The slopes of the 
edges of this region are related to the critical exponents according to Campi [|14[ . 



Clearly there are differences between percolation theory and a Gibbs model, but the dif- 
ferences are not as large as originally suggested by early computations. Indeed the reduced 
variance might not reliably distingush percolation from a simple Gibbs model. A different 
method is needed to distingush these models. However, the idea of excluding the largest 
cluster from the ensemble averages is a standard procedure in percolation theory flT5| , and 
this new technique for doing that analytically in the Gibbs models shows a particular advan- 
tage of these models over percolation models, which we hope will encourage further interest 
in them. 

This work supported in part by the National Science Foundation Grant No. NSFPHY 
92-12016. 
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FIGURES 

FIG. 1. Expected reduced multiplicity (72} (m) (a) and largest cluster size k max /A (b) vs. 
(m — 1)/A for r = 1.0,2.5 and the percolation model at A = 125. 

FIG. 2. Campi probability contour plot for Z = 79, r = 2.5 at the critical point x = x c . 
The axes are logarithmic, with the largest cluster size on the y-axis, and the ratio of the reduced 
moments M2/M1 on the x-axis. The central rings are higher in probability than the outer rings. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9508038vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/nucl-th/9508038vl 



